Optimal dynamics for quantum-state and entanglement transfer through 

homogeneous quantum wires 
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It is shown that effective quantum-state and entanglement transfer can be obtained by inducing a 
coherent dynamics in quantum wires with homogeneous intrawire interactions. This goal is accom- 
plished by tuning the coupling between the wire endpoints and the two qubits there attached, to an 
optimal value. A general procedure to determine such value is devised, and scaling laws between the 
optimal coupling and the length of the wire are found. The procedure is implemented in the case 
of a wire consisting of a spin-^ XY chain: results for the time dependence of the quantities which 
characterize quantum-state and entanglement transfer are found of extremely good quality and al- 
most independent of the wire length. The present approach does not require ad hoc engineering of 
the intrawire interactions nor a specific initial pulse shaping, and can be applied to a vast class of 
quantum channels. 
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One of the most commonly requested conditions in 
quantum communication and computation protocols is 
that two distant parties, typically Alice and Bob, share 
a couple of entangled qubits. When the physical objects 
encoding the qubits can travel, as in the case of optical 
photons, the above goal can be accomplished by creat- 
ing the entangled couple in a limited region of space and 
then letting the qubits fly where necessary. On the other 
hand, when qubits are realized via intrinsically localized 
physical objects, as in the case of 5 = ^ spins or atomic 
systems, a different strategy must be adopted (see for in- 
stance Ref. J, and references therein). One such strategy 
is the following: first, two neighboring qubits (A and A') 
are prepared in an entangled state, by means of a short- 
range interaction; then, the mixed state of one of the two 
qubits (say A) is transferred to a third distant qubit via 
a quantum channel. If state-transfer is perfect the proce- 
dure results in a pair of distant entangled qubits A' and 
B, as requested. 

Aim of this paper is to set a general framework where 
such strategy can be successful. In particular, {%) we 
define a procedure for controlling such dynamics, and 
hence the quality of the transmission process, by spe- 
cific operational settings; (n) we show that the quality 
of the quantum-state and entanglement transfer is not 
substantially affected by the length of the wire; {Hi) we 
apply the procedure to the spin-^ XY chain and show 
that high-quality quantum-state and entanglement trans- 
fer are obtained. 

Let us first recall that for the strategy depicted above 
to make sense, one has to equip oneself with a quantum 
channel capable of transferring mixed states. How to 
obtain such a channel is the problem to which many au- 
thors have proposed different solutions (Tj-Q ? some based 
on the idea of engineering the channel itself, by the spe- 
cific design of its internal interactions, others on that 
of intervening on the initialization process, by preparing 
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FIG. 1: The endpoints of a quantum wire V are coupled to 
the qubits A and B, via a switchable interaction j; A can be 
entangled with an external qubit A^ 



the wire in a configuration found to serve the purpose. 
In both cases, a severe external action on the physical 
system is required. 

Here a different point of view is adopted: the pur- 
pose is that of devising conditions for an optimal dy- 
namics to occur, where by 'optimal dynamics' we mean 
a time evolution of the state of the channel such that the 
mixed state of A is transferred with high fidelity through 
the channel to a distant party B. The main guideline is 
the fact that excitations characterized by a linear disper- 
sion play a crucial role in determining effective transfer 
of quantum states. Since this idea has been put for- 
ward [l, m, Q, it had to face the unfortunate evidence 
that a quantum channel with a linear dispersion relation 
over the whole Brillouin zone is not only difficult to de- 
sign [4], but perhaps just a chimera to realize, so far. On 
the other hand, most physical systems are characterized 
by excitations whose dispersion relation has zones of lin- 
earity, typically near inflection points: if we were able to 
induce a dynamical evolution of the channel essentially 
ruled by excitations belonging to such zones, a coher- 
ent propagation should result, and an effective transfer 
of quantum states consequently obtained. Let us hence 
focus on how to induce such dynamical evolution. 
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We consider (see Fig. [T]) two separated parties, Alice 
and Bob, each owning a qubit, A and B, respectively. 
Each qubit may be subjected to a local interaction, h^l-L^ 
and h^l-L^^ respectively, with possibly tunable. The 
quantum channel is realized by an extended physical sys- 
tem r, hereafter called quantum wire, which is here as- 
sumed to be made of N interacting particles on a dis- 
crete lattice. The internal dynamics of the wire is ruled 
by the Hamiltonian H^. The endpoints of the wire can 
be put in contact with the qubit A (B) via a switchable 
interaction jT/^r (j^sr)- The overall system has mirror 
symmetry and the Hamiltonian which rules its dynamics 
is H = + j + '^Br ) + a + '^b ) • Before the 
process starts A, B and T do not interact with each other; 
the wire is in its ground state l^^r)^ while Alice and Bob 
prepare their qubits in the initial states Pj^{0) and Pj^{0), 
respectively. At time t = the interaction between each 
qubit and the wire is switched on. 

We assume that: (i) H can be written (either exactly 
or via a reasonable approximation) as a quadratic form 
of N -\- 2 local, either fermionic or bosonic, operators 
{r]i, 7^|}, where the index i = 0, . . . , N-\-l labels the sites 
of the lattice; (ii) the bulk of the wire is translation sym- 
metric. Condition (i) guarantees the existence of a uni- 
tary transformation {rji^rjl} {Vk^vl}^ which diagonal- 
izes the total Hamiltonian, H = "^^^kvlvk- Condition 
{ii) entails that the above transformation is close to a 
Fourier transform when i corresponds to sites in the bulk 
of the wire, so that k can be approximately considered 
as a quasi-momentum. For the sake of clarity, let us 
consider the qubits A and B as initially prepared in the 
pure states \a) and respectively; the overall system 
at time t = is then described by |^o) = \(^)\^r)\P)^ 
which generally is a non-equilibrium state. Its dynamics 
is ruled by quasi-particle excitations with density in k- 
space n{k) = (^o|^i^/c|^o), which, in their turn, evolve 
according to H. If the dispersion relation were linear, 
ujk — OjQ -\- vk, the bulk dynamics would consist in a coher- 
ent wavepacket traveling from A to B (and vice versa) at 
velocity v along the wire. Under our further assumption 
of mirror symmetry, such wavepacket would perfectly re- 
build the initial state of A on the qubit B, after a time 
tc^N/v. 

However, dispersion relations of real interacting sys- 
tems are typically non linear, and display a much more 
complex dependence on the wavevector k (see for instance 
Fig. [2]). On the other hand, if we cannot deal with a most 
unusual linear dispersion, yet we can act on the distribu- 
tion n{k) by varying the coupling between the qubits and 
the wire, or other parameters entering the Hamiltonian 
of the qubits, or their initial states \a) and \f3). 

Suppose then that we can adjust the distribution n{k) 
to make it peaked around an inflection point ko of the 
dispersion relation, thus setting the most relevant exci- 
tations for the dynamics, as those whose energies can be 
written as 

(^k ^ ^/co + v{k-ko) + y (/c-^o)^ , (1) 
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FIG. 2: Dispersion relation of the spin-^ XY chain, for dif- 
ferent interaction parameters. The bell-shaped curve edging 
the shaded area is the optimal density of excitation vs. cj, as 
obtained for 7 — 70, at — 0.5, and h — 0.5. 

where the main dispersive term is retained. Further sup- 
pose that we can set the width of the distribution n(/c): 
be 1/(4(7^) its variance. It can be shown [10] that, under 
the above conditions, the dynamics of the wire is still 
described by a traveling spacial wavepacket, centered at 
x(t) = (v -h acr~^) with a variance that increases from 
the initial value (7^(0) = as (7^(1) = + aH^I(2(j'^) . 

We hence state that the optimal dynamics is obtained 
by choosing so as to minimize cr(t) at the arrival time 
t N/v, leading to 

<Topt= (|a|iV/t;)'/' . (2) 

As stated above, the form of the dispersion relation 
is essentially designed by the wire Hamiltonian, the po- 
sition of the peak of n{k) can be set by acting on h^, 
while (J is expected to depend on the coupling between 
the qubits and the wire, whose energy scale is j. 

Based on the above reasoning, we define the follow- 
ing procedure for inducing an optimal dynamics: Fix the 
Hamiltonian parameters of the wire so as to grant the 
existence of a wide enough interval in k with almost lin- 
ear dispersion as described by Eq.([T]); tune so as 
to peak n{k) at an inflection point of uok; then set a to 
its optimal value Eq. (|2]) by varying j. The value of the 
coupling corresponding to (Topt will be hereafter referred 
to as the optimal coupling, jopt- It is worth mention- 
ing that the optimization of a can also be obtained by 
creating a finite-size excitation in space as the starting 
state [3, [5I, [m, [l3|; however, this strategy results in a 
much less versatile procedure. The weak dependence on 
N embodied in Eq. ([2]) is expected to give small losses for 
transmissions over long wires. In fact, the actual time t^ 
needed for the packet to travel along the wire, implicitly 
defined by x{tN) = is found to depend on N accord- 
ing to vtN — N — (Topt . It is of absolute relevance 
that cr(tAr) — Y^3/2 (Jopt , i.e., the final optimal packet is 
just about 22 % wider than at start, irrespectively of the 
wire length. 

Let us now describe a specific implementation of the 
procedure described above. We consider a system de- 
fined on a one-dimensional discrete lattice of N-\-2 sites. 
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labeled by the index i = 0, 1, . . . , The qubits 

A and B sit at sites and A/'+l, respectively. The wire 
is physically realized by a chain of N interacting S =^ 
spins, taking up sites from 1 to TV. Neighboring spins in- 
teract via a Heisenberg Hamiltonian of the XY type, and 
are possibly subjected to an external magnetic field. The 
exchange interaction and the magnetic field are assumed 
homogeneous along the wire. The total Hamiltonian is 

N-l 

n= - ^(l+7)5f5f+i + {l-j)SfSf^,-h J2 St 

i=l i=l« 

-3 E [(l+7.)5f + (l-705f 

-/iJ5o' + 5^+i), (3) 

where the exchange energy for the intrawire interaction 
has been taken as the reference energy scale and hence set 
to unity. Mirror symmetry implies =70- An isotropic 
exchange interaction in Eq. O 7 = 70 = 0, defines the so 
called XX model. A Jordan- Wigner transformation, 
casts the Hamiltonian (j3j) into a quadratic form of 
interacting fermionic operators, each defined on 

a site of the lattice. A further Bogolubov transforma- 
tion diagonalizes the Hamiltonian, 1-L = J^/c ^kvlvk + ^0, 
where £^0 is the ground state energy, whereas 7^^ {r]k) are 
fermionic operators which create (annihilate) excitations 
of energy uJk [14]. 

The optimization procedure described above begins 
with the analysis of the dispersion relation in the infinite 
chain limit, cok = [{h— cos k)'^ + 7^ sin^ ^]^^^, displayed in 
Fig. [2j One can easily spot the existence of more or 
less wide regions of linearity in the neighborhood of the 
infiection point (s) /cq, where Eq. ([T]) holds. By tuning 
— UJkoi the peak of n{k) is made to sit at k^. The op- 
timal coupling jopt is then numerically determined so as 
to fulfill Eq. (|2]). In fact, the dispersion relation strongly 
depends on the parameters of the wire Hamiltonian, 7 
and h: in particular, the region of linear dispersion sen- 
sibly shrinks as the anisotropy 7 increases, which might 
make the wire to be useless; however, the linear region 
can be extended again by increasing the field h\ there- 
fore, one can act on the latter parameter so as to fulfil 
the conditions for optimal dynamics. For example, in 
the extreme case of the Ising chain (7 = 1) for /i = the 
dispersion relation becomes fiat and does not allow for 
propagation; however, a wide linear region can yet be 
obtained by applying a finite h onV. 

Let us now give a fiash upon the resulting dynamics of 
the overall system. To this aim, we set 7 = 70 = {XX 
model) and /i = 0, as this choice makes more analytical 
expressions available, which in turn allows for a more de- 
tailed analysis. The infiection point /cq = 7r/2 corresponds 
to uOkQ = and we hence set h^=0. 

In Fig. [3] we show the time evolution of the magne- 
tization parallel to the quantization axis all along the 
wire, (5'f(t)),z = 0, 1, . . . , TV, A/'+l, for N = bO and ini- 
tial states of the qubits |a) = = The upper panel 
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FIG. 3: Time evolution of the on-site ^-magnetization for 
N =50 and an initial state |t)|^r)|t); the qubit-wire coupling 
is J = 1 (top) and j = jopt = 0.58 (bottom). 



corresponds to a generic value of the coupling, j = 1. The 
lower panel is for j=jopt = 0.58: this value is deter- 
mined exactly by exploiting analytical expressions hold- 
ing for the XX model [l5[. The difference from the pre- 
vious panel is striking: indeed, the dynamics of the wire 
is essentially ruled by the nondispersive propagation of 
wavepackets, and results in {S^{t)) being, at the arrival 
time t an almost perfect reproduction of the initial 

magnetization {S^{0)) of the qubit A. These results do 
also confirm that our prescription for the above coherent 
propagation to occur, i.e., for the determination of jopt, 
is correct. 

But does this peculiar dynamical evolution also affect 
the quality of quantum-state transfer? To answer this 
question, we now deal with the time evolution of quan- 
tities which are used to monitor such quality. In par- 
ticular, we analyze the quantum-state transfer process 
in terms of the fidelity between an initial state \a) of A 
and the evolved state pj^{t) of B, J^j^^{t) = {a\p^{t)\a), 
where initially Pb(0) = lt)(t|- As for the entanglement, 
we refer to the time evolution of the concurrence between 
A' and B, C^j^,{t) = C{p^^,{t)) [16]. In the case of the 
XX model we can also evaluate the minimum fidelity 
over all possibile hereafter indicated by J~'^^^{t), 
which allows us to put forward further conclusions about 
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FIG. 4: Time evolution of the minimum fidelity J^™g^(t) vs j, 
for aA/' = 50, 7 = 7o=7Ar = {XX chain) , and h — h^ =0. 



the quality of the entanglement transmission from AA' 
to BA^ Indeed, this can be done using the following 
lower bound [13, Ell for the fidelity of entanglement: 

•^(I^aa')(^aa'Uba'W) > l-^re^W-i The entangled 
state 1^^^/) can be hereafter taken as any of the Bell 
states, since a local operation on A' does not change the 
dynamics of the concurrence. 

The above quantities essentially depend on the time 
evolution of the magnetizations of B, which we have de- 
termined as follows: The components have been 

written in terms of the fermionic operators {^/c^^}; 
their time-dependence, which simply follows from 1-L = 
^k^kVkVki has been made explicit, so as to obtain, by 
transforming back via inverse Bogolubov and Jordan- 
Wigner transformation, the Heisenberg representation of 
the spin operators S^j^-^{t). Their expectation values, 
i.e., the required magnetizations, have been finally de- 
rived by expanding S^j^-^{t) in terms of the spin op- 
erators of the extremal qubits and of fermionic oper- 
ators relative to the wire, i.e., defined by the diago- 
nalization of the sole This approach allows us to 

devise a numerical procedure for determining, for any 
model belonging to the class described by the Hamil- 
tonian (|3]), the Kraus operators |17|] M^(t), in terms of 
which we get = E^^MWpA(0)Mt(t) and p^^, (t) = 
T.^, [M^,{t) l]pAA,(0)[Mt(t) l] , i.e., the necessary 
tools for our analysis. 

We first consider the fidelity J^^^ (t) as defined above. 
In the case of the XX chain, we put ourselves in the worst 
possible case and evaluate J^'^^it)'- a high value of such 
quantity ensures a very good transfer of any initial state, 
modulo a local operation. In Fig. |4]we show ^^^"(^) for 
the same model as in Fig. [3] as a function of the coupling 
j. A minimum fidelity which is just slightly below unity is 
observed for that same value of the coupling, determined 
by our procedure, which gives rise to the dynamics shown 
in the lower panel of Fig. O i.e., j = jopt — 0.58. 

We underline that such a high value of the minimum fi- 
delity, which is confirmed also for as large as 500, does 
also imply an optimal transmission of entanglement. It 
is worth noticing that the peak of J^^^it) occurs simul- 
taneously with that of {S^{t))^ i.e., at t ~ A', a feature 




FIG. 5: Time evolution of Tab (upper panel) and C^^, (t) 
(lower panel), for different A, 7 = 70 = 0.5, h = 0.5, /Iq = 0.85, 
and jopt= 0.49, 0.39, 0.34 for A = 50, 250, 500, respectively. 
The second signal observed in both panels is the first echo 
for A = 50. 



that can be analytically proven in the XX case. 

The consistence of the above picture allows us to move 
forward to the analysis of the entanglement transfer, 
and also in the more general case of the XY model. 
We set 7 = 70 = 0.5 and h = 0.5 and apply our proce- 
dure: We first locate the inflection point of the disper- 
sion relation in the thermodynamic limit at /cq — 1.795 
and tune h^c^O.Sb. This latter value slightly differs 
from cj/cQ as determined in the thermodynamic limit, due 
to finite-size and boundary effects. Finally, Eq. (|2]) al- 
lows us to numerically determine jopt = 0.49, 0.39, 0.34 for 
A" = 50, 250, 500, respectively. 

Despite the dynamics in the XY case being compli- 
cated by the existence of several energy scales, we man- 
age to get a very good state transfer also for large A", 
as testified by the fidelity, averaged over all possible ini- 
tial states of A, /"^^ shown in the upper panel of FiglS] 
The oscillations observed are due to the precession in- 
duced by the magnetic field acting on the qubit B. 
Such precession causes only a phase change and does not 
affect the entanglement transfer, as shown below. The 
entanglement between B and A' is monitored by the time 
evolution of the concurrence C^^,{t)^ shown in the lower 
panel of Fig. [5] for different values of A^: For each A', 
this is characterized by an extremely well defined peak, 
of height > 0.8 even for chains as long as 500 sites. This 
peak occurs simultaneously with the maximum of 
and has a finite but small width. 

The second signal observed in both panels is the first 
echo for A" = 50, which is seen to be still considerably 
intense and well localized: this testifies of a long-lasting 
quasi non dispersive dynamics of the wire. Other echoes 
occur at times which are too large for being shown in the 
figure, but their structure confirm this statement. 

In conclusion, based on a general picture of the dynam- 
ical evolution of quantum wavepackets, we have devised 
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a procedure for inducing an optimal dynamics in a quan- 
tum wire. The procedure leads to the determination of 
an optimal coupling between the wire and the two qubits 
attached at its endpoints, meanwhile giving indications 
about the best setting of other, possibly tunable, param- 
eters of the system. For the procedure to apply, few very 
simple conditions must be fulfilled, and there is no need 
for a specific design neither of the wire, nor of its initial 
state. By implementing our approach to the spin-^ XY 
chain, we have obtained extremely good quantum-state 
and entanglement transfer. The time scales over which 
such transfer is found to occur is considerably shorter 



as compared with previous results concerningquantum- 
state transmission over Heisenberg models [l|, [HI, 
21]. Moreover, the quality of the state and entanglement 
transfer that we obtain very weakly deteriorate as the 
length of the wire increases. 
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